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We consider the dependence of the MgE>2 superconduct- 
ing critical temperature on the pressure. Our model exploits 
the influence of the large polarons on the band structure of 
the layered MgB2 superconductor. Namely, the hole Pekar- 
Frohlich polarons form quasi two-dimensional potential wells 
in the boron plane which shift the positions of the a- and 
7r-bands. This energy shift depends on the pressure and the 
Cooper pairing of the correlated er-electrons happens inside 
polaron wells. The results obtained are as follows: dT c /dp ~ 
-a(5.2 ± 0.9) K/GPa or dT c /dp ~ -a(6.9 ± 1.1) K/GPa for 
a different choice of the Griineisen parameter. Being com- 
pared with known experimental data they give us a resonable 
interval for the value of the Frohlich electron-phonon coupling 
constant: a ~ 0.15 — 0.45. 

PACS numbers: 71.38.Fp, 74.25.Kc, 74.62.Fj, 74.70. Ad 



I. INTRODUCTION 

Soon after the discovery of the superconductivity in 
MgB2 the pressure influence on the critical tempera- 
ture T c has been experimentally studied. The baric 
derivatives dT c /dp were extracted from the obtained 
data. At the low-pressure (p < 2 GPa) the experi- 
mental results are as follows: dT c /dp = —1.2 K/GPa 
in Refs. U and -2 K/GPa in Refs. |,||. Exper- 
iments at high-pressure lead at the same sign of the 
baric derivatives but demonstrate lower values. Namely, 
dT c /dp = (-0.8,-1.07, -l.ll)K/GPa in Refs. §011 re " 
spectively. Many conventional low-T c superconductors 
also reveal the same sign, dT c /dp < 0, but their absolute 
values are much lower. 

Borons in the crystal MgB2 are packed in honeycomb 
layers alternating with hexagonal layers of magnesium 
ions. The in-plane distance, b = 1.78 A, between the 
boron cations is less than those between magnesium an- 
ions, a = 3.08 A. The space separating the boron planes 
is of the size c = 3.52 A. The magnesium ions are posi- 
tioned above the centers of hexagons formed by boron 
sites and donate their electrons to the boron planes: 
Mg 2+ [B~ (p 2 )] 2 - These p-electrons form cr- and 7r-bands 
and the charged Mg-layers lift the three-fold degeneracy 
between a- and 7r-electrons and shift the electronic en- 
ergy bands so that 7r-band becomes lower than the o- 
band. It occurs that the latter crosses the Fermi-level 



providing the light and heavy hole formation. The elec- 
tronic structure is formed by the narrow energy-band of 
the two-fold degenerate a-electrons and the wide-band of 
the 7r-electrons. 

As to the phonon modes in MgB2 they have a sharp 
cut-off at about 100 meV. Below this energy it. is sup- 
posed the existence of different phonon modesQ. Note 
that the phonon energies mentioned in Ref. [?] could be 
overestimated because the approximation used to derive 
them is valid only for cubic latticesEm. The anharmonic 
in-plane breathing mode Ei g (74.5 me VX. strongly inter- 
acts with the narrow band <7-electronsEl. Due to this 
fact its energy is lower than that of the out-of-plane 
tilting boron mode B\ s (87.1 meV). Besides, there 
are lojaij-eafirgy acoustic modes and indications are also 
foundEftj of the existence of phonon modes with ex- 
tremely low energies. 

The most of the theoretical speculations on the T c 
pressure dependence in MgB2 are based on the semi- 
phenomenological approaches and do not involve the mi- 
croscopic arguments. In the scope of the conventional 
BCS theory dT c /dp is governed by the competition be- 
tween a decrease of the density of electronic states at the 
Fermi level and an increase of the Debye frequency in 
the BCS formula for T c . The weak electron-phonon cou- 
pling constanto in MgB2 and the absence of the peak 
in the electronic density of states at the vicinity of the 
Fermi-energyll3 cannot provide the BCS scenario. 

On the other hand, the spectroscopic measurements 
of the superconducting gaptjlla,— ..the specific heat 
behaviorE3E3, the low isotope effecttilo (under the sub- 
stitution of_both B and Mg isotopes) and the pres- 
sure effectsEJ evidence in favor of the complex mech- 
anism of the superconductivity in MgB2. | —.The calcu- 
lated spectral functioned, and the analysisEj of the re- 
flectance measurementsEj show the possibilities of dif- 
ferent superconducting mechanisms beyond the conven- 
tional electron-phonon BCS pairing. 

The positive baric derivative of T c was foundElEU in the 
microscopic model of the hole dressed superconductivity 
what opposes the experimental data. 

Among optical modes in MgB2 the Ei g is the only Ra- 
man active mode. r-The high-pressure (up to 15 GPa) 
Raman experimenter has revealed a large pressure shift 
of the Ei g phonon mode. The Eig mode does not man- 
ifest any softening at T c (see Ref. p9| ) therefore does not 
provide the superconducting electron pairing. Neverthe- 
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less, it is responsible for the T c behavior with pressure as 
we will see later. It should be pointed out that baric ex- 
periments have been completed both for pressed powder 
and single crystals_No structural transitions were found 
in experiments0c3'E2l up to the pressure 40 GPa. 

The goal of the present paper is to calculate the baric 
derivative dT c /dp. We suppose that the large polarons 
are formed in the boron-planes due to the er-electrons 
interaction with the optical E^-phononsQ. Because of 
the polaron formation the <r-band becomes lower while 
the 7r-electrons are not influenced. As the events take 
place in the layered medium we consider the quasi-two- 
dimensional polaron. The polaron anisotropy parameter 
and, correspondingly, the er-band energy shift depend on 
the lattice constants influenced by the pressure. Roughly 
speaking, the Cooper pairs of the correlated cr-electrons 
in the MgB2 superconducting condensate are plunged 
into the pancake polaron. To treat the Cooper paring 
we accept the kinematic mechanismEirE3 in MgB2. 



II. FORMULATION OF THE MODEL 

In order to describe the superconducting state of 
MgB 2 , we adopt the model with the strongly interacting 
cr-electrons (with a narrow bandwidth 2w\) overlapping 
with the noncorrelated 7r-electrons (with a wide band- 
width 2u>2) of boron ions. The important parameter here 
is the distance r between the centers of a- and 7r-bands 
influenced by the pressure. 

The positive-Hall current and the thermoelectric power 
measurements^ indicate the hole conductivity in MgB2. 
The polaron radius is given by the expression 
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where m and m e are the hole and the vacuum electron 
masses, respectively, and fuo is the low-frequency opti- 
cal phonon energy responsible for the polaron formation. 
Inserting the value fttu ~ 74.5-^7161^ for the _E 2 g-phonons 
and the light hole band mass£3 m — 0.25m e into Eq. (|l|) 
the polaron radius can be estimated as R ~ 20 A which 
is much larger than the lattice constant. 

Studying the polaron sector of the system we apply the 
conventional Pekar-Frohlich Hamiltonian generalized for 
the consideration of the peculiarities of MgB2. To take 
into account the fact that the cr-electrons are located in 
narrow layers inhabited by boron atoms we introduce the 
in-plane mass m which is the light hole band mass and 
the large mass M in the orthogonal (z) direction. The 
ratio m/M of these masses is related to the probability 
of a carrier to tunnel from one layer to another. 

To take into account the different polarizability of the 
medium in the rcy-plane and in the z-direction we sug- 
gest to simulate it by using the anisotropic Columb po- 
tential for the interaction of the hole with the induced 
polarization field. That is, we use the potential of the 



form V(r) cx 1/ ' y/ x 2 + y 2 + (z/7) 2 with 7 < 1 being a 
phenomenological parameter to describe the deformation 
of the Coulomb potential along the z-axis (that is, the 
oblation of the equipotential surface) . In the momentum 

space the potential is of the form V(k) cx l/(fc| + j 2 k 2 ). 

In a sense we study the electronic sector as how all 
events happen to be located inside the large polarons. 
The large polaron being formed we obtain the lowering 
of the energy of carriers, that is the renormalization of 
the parameter r which becomes dependent on the applied 
pressure. This influences drastically the electronic struc- 
ture of MgB2-like systems, especially the relative position 
of a- and 7r-bands. 



A. Polaronic Sector 

Thus, we consider an anisotropic polaron whose motion 
is confined to the xy-plane while the motion in the per- 
pendicular z-direction (along the crystallographic c-axis 
of MgB2) is restricted. The Pekar-Frohlich Hamiltonian 
reads as follows: 
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where fl is the crystal volume, p\\,p z are the hole mo- 
menta in the xy-pl&ne and in the orthogonal z-direction, 
respectively, and uj is the LO-phonons frequency. Oper- 
ators at(ctg) are the creation (annihilation) operators of 

the phonons with a momentum k. The electron-phonon 
interaction term is specified by Vr; its squared modu- 
lus is the Fourier transform of the Coulomb potential 
V[f). In contrast with the conventional definition we in- 
troduce an additional anisotropy of the Coulomb poten- 
tial to describe the different polarizability of the medium 
in various directions (in the isotropic boron plane and 
perpendicular to it): 
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Here a is the conventional Frohlich coupling constant of 
the electron-phonon interaction. nr-i|— 1 

The MgB2 baric compression is anisotropicatflO. Ac- 
cording to Ref. |28| the compressibility along c-axis al- 
most twice larger the plane compressibility. Under the 
hydrostatic pressure the initial compressiop-|along the c 
axis is larger than along the boron planec3. We esti- 
mate the parameter 7 as the ratio of the standard lattice 
constants: 7 ~ c/2a (here c/2 is the distance between 
charged boron and magnesium planes). When the pres- 
sure increases the distance between the Mg and the B 
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planes decreases which is described by the decreasing of 
the Coulomb potential anisotropy 7 in our model. The 
polaron self-energy AE can be found within the second 
order of the perturbation theory: 



AE = - 
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The integration in Eq 
the expression: 



being performed we arrive at 
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Here the parameter T = 7^/Af/m. This formula repro- 
duces two well-known limiting cases. When M = m and 
7 = 1, then r = 1 and we obtain A_B = —ahuj, that 
is the conventional result for the three-dimensional po- 
laron. When 7=1 and M 3> m, we arrive at the result 
AE = —jafiio which is vaJid .for the two-dimensional po- 
laron confined to a planeOEa. If 7 is finite and M 3> m, 
then f>l and Eq. (|J) leads at the following expression: 
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It is the MgB2 case indeed on the reason that the mass 
ratio is inversely proportional to the hopping integrals 
ratio: M/m oc t/t z oc 10 (see Ref. ^). For the nu- 
merical estimates we will take the experimental value for 
the optical phonon frequency to. The energy shift AE is 
negative what means the total lowering of the band min- 
imum of the hole carriers due to the polaron effect. Due 
to AE the energy distance r between the a- and 7r-bands 
is shifted by a frequency dependent contribution which 
is essential for the electronic sector. 



B. Electronic Sector 

We start here from the high-temperature paramagnetic 
phase of MgB2 system and derive the superconducting 
critical temperature from the condition of instability of 
the normal state of correlated cr-electrons with a tem- 
perature decrease. Namely, the T c is governed by the 
solution of the Bethe-Salpeter equation for a vertex Y p 
in the Cooper channel in the reference frame of the elec- 
tron pair: 
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where G°J (q) = l/(- 



■IU) T , 



■ £ q ) is the normal state 



Green's function for a strongly correlated cr-electron with 
a spin orientation s and an.energy dispersion £ g = ft q — (i 
with the correlation factored 
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The Matzubara frequencies are given by uj n — (2n+l)nT 
inEq. (|). 

The expression for the chemical potential 
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follows from the equation n a + n w = 2 for the total elec- 
tron density per boron at the assumption of a complete 
ionicity of Mg ions in Mg 2+ [B~ {p n p +n ")] 2 system. In 
the integral kernel of Eq. (Q) near the Y-A line of the 
Brillouin zone the Coulomb vertex V (p — q) can be fac- 
torized as 



V(p-q) = 2/3t p t q , 



(10) 



where the parameter (3 — V/6t 2 labels an effect of the 
Coulomb repulsion for the nearest cr-electrons and the 



energy dispersion is t p = 3i 1 1 
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p ) /12] near the 



Y-A line. The kernel of the integral equation (^) does 
not contain another kinematic verticesEJ which could be 
essential for superconducting condensates at moderate 
densities of carriers and/or for specific symmetries of su- 
perconducting order parameters. 

Summation over the Matzubara frequencies being per- 
formed, Eq. (0) takes the form: 
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The superconducting coupling constant can be written 
as 
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so the chemical potential is restricted by the inequality 
< £t — / / P which generates constraints for the elec- 
tron structure parameters and r. The larger is the 
electron-electron Coulomb repulsion V ~ /3 the more nar- 
row becomes the superconducting region and the lower 
is T c . 

Hereafter we will neglect the inter-electron Coulomb 
repulsion (/3 = 0). Then the superconducting critical 
temperature satisfies the following equation: 



(13) 



where the £(wi) and £(—wi) are the energy dispersion 
values of cr-electrons at the top and the bottom of their 
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energy band, respectively, and p(£) is the electronic den- 
sity of states. Eq. ( |l3| ) for the superconducting critical 
temperature coincide with one derived in Refs. |3^j3^ in 
a different way. 

Note that the kinematic superconducting mechanism 
was also applied to MgB2 in Ref. fy] but with a non- 
physical negligence of the 7r-electrons role. In our ap- 
proach the great importance has the characteristic en- 
ergy difference r between of a- and 7r-bands, shifted by 
AE due to the interaction between the light holes in the 
cx-band and the quasi-2D E 2g phonons: 



r + AE. 



(14) 



Due to the chemical bonds nature the hydrostatic 
pressure decreases the inter-plane distance more read- 
ily than the in-plane boron-boron distance, so that 
dw 2 /dp 3> dwi/dp (c.f. Ref. [l2]) and the latter deriva- 
tive can be neglected in our estimates. The calculations 
of the T c (r) demonstrated that r > —Wi/4 for MgB2. 
As the next step we differentiate the integral equation 
( |l3| ) with respect to the pressure and take into account 
Eq. (|l4|). Assuming the rectangular density of states 
p = 9{w\ — ^ 2 )/2wi we obtain after the subsequent inte- 
gration the following expression: 
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where 70 = InC with C = 0.577 being the Euler's con- 
stant. Under the natural assumption T c -C w\ one can 
keep only the first term in the brackets in Eq. (15). Then 
we arrive at the result: 
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Taking into account that the Frohlich coupling con- 
stant a oc 1/V^ we fi n( i then from Eq. (^): 
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dp 



dp 
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The parameter 7 = c/2a evidently decreases with the 
pressure. As to the quasi-2£> frequency lo we use the esti- 
mate of a quasi-harmonic phonon via the mode Griineisen 
parameter d\nuj/dp = G/Bq. The quantity Bq is said 
to be the bulk modulus. It follows then that both w 
and \/^/7 increase with the pressure. Tliesp-findings are 
compatible with an experimental studycSo of Raman 
spectra and lattice parameters of MgB2 under pressuxe 
at room temperature and with theoretical estimations^ 
of the strong influence of the pressure on E% g mode. Fi- 
nally we conclude that dAE /dp < and dT c /dp < 0. 

In an explicit form the baric derivative for the super- 
conducting critical temperature can be written as follows: 
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III. NUMERICAL RESULTS AND 
CONCLUSIONS 



For numerical calculations we estimate the polaron 
anisotropy parameteii-as 7 = c/2a = 0.57. The known 
value of its derivativeE3 is d\n~/ 2 /dp — 2d\n (c/2a) /do ~ 
—2.4 ■ 10~ 3 GPa^ 1 . The bulk modulus is measuredM to 
be B ~ 114 GPa. The Griineisen parameter reported 
in Ref. [2^ equals G = 2.9 ± 0.3 for the measured Raman 
active i?2g-phonon mode with the energy huj = 76.7 meV 
which we will use in Eq. (fL8|). In particular, the polaron 
self-energy AE = —211.4a meV. The realistic value for 
the MgB2 energy bandwidths ratio of n- and tr-electrons 
is w 2 /wi ~ 18/9 = 2 (see Ref. 

Putting these magnitudes in Eq. (|l8[) we obtain the 
estimate for the derivative of T c with respect to pressure: 



dT c 
dp 



-a(5.2±0.9) 



K 
GPa' 
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The uncertainty of the result comes from the experimen- 
tal deviation errors in the Griineisen pata-meter. The 
comparison with the low pressure resultlM dT c /dp = 
-1.2 K/GPa leads at the interval a = 0.20 - 0.28 for 
the Frohlich electron-phonon-pCoupling constant. The 
comparison with the resultldoa —2 K/GPa leads sub- 
sequently at the estimate a = 0.33 — 0.47. 

These results are shown in Fig. 0. The solid line 
presents the dependence of the baric derivative dT c / dp on 
the electron-phonon coupling constant a. The thin solid 
lines set the experimentally measured values for dT c / dp 
and the dashed lines show error bars for a at a taken 
value of the baric derivative. 

We have to note that the numerical results can be 
changed if one scales the frequency shift of the in-plane 
mode with the variation of the interatomic bond distance 
or lattice parametercZI. Then, the Griineisen parameter 
takes even the larger valuec3 G = 3.9 ± 0.4 and our 
Eq. (|l9| ) takes the form: 
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Subsequently, the estimates for the electron-phonon cou- 
pling constant will be different as well: a = 0.15 — 0.21 
for dT r Jdp = -1.2 K/GPa and a = 0.25 - 0.34 for 
dTc/dp = -2 K/GPa. 

Combining the numerical results for a we may con- 
clude that the electron-phonon coupling constant in 
MgB2 is in the range a = 0.15 — 0.45 which seems to be 
quite reasonable interval for the assumed weak-coupling 
regime of a large hole polaron. 

To resume, we proposed a model in which large 
anisotropic polarons play an important role decreasing 
the energy distance r between the a- and 7r-bands. The 
polaron anisotropy is governed by the introduced param- 
eter 7 which depends on the geometry of the MgB2 sys- 
tem influenced by the pressure. Besides, the phonon fre- 
quency increases with the pressure which results also in a 
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decrease of the polaron radius. Thus, the superconduct- 
ing properties of the MgB2 system are influenced both 
by the geometry of the crystal and by the polaron well 
depth and size. 

The superconducting instability is driven by the non- 
phonon kinematic mechanism with T c depending on the 
energy difference r between the a- and 7r-bands. The 
quantity r incorporates all electron-phonon effects in our 
model. As the result the baric derivative dT c /dp is calcu- 
lated as a function of the Frohlich coupling constant a. 
At realistic values of a our calculations agree with the ex- 
periments. It follows from the arguments discussed in the 
present paper that the pressure measurements provide us 
with valuable tests to establish the acceptable model for 
MgB2, to understand the mechanism of the supercon- 
ductivity in this material and to estimate its Frohlich 
electron-phonon coupling constant. 
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FIG. 1. Dependence of the baric derivative of the supercon- 
ducting critical temperature on the Frohlich electron-phonon 
coupling constant (solid line). The dashed lines show er- 
ror bars iof a due to the uncertainty in the Gruneispp, 
parameteru. The thin solid lines show experimental valuesBu 
for dTJdp. 



ACKNOWLEDGMENTS 

We are grateful to K. Kikoin for stimulating discus- 
sions. M.A.S. thanks Universiteit Antwerpen — TFVS 
for the hospitality during his visit to Belgium. The work 
is supported by GOA BOF UA 2000, IUAP and FWO-V 
project No. G.0274.01. 



General and Inorganic Chemistry of the Russian Academy 
of Sciences, Leninskii prospekt 31, 117 907, Moscow, Russia 
* Also at: Universiteit Antwerpen (RUCA), Groenenborger- 
laan 171, B-2020 Antwerpen, Belgie and Technische Uni- 
versiteit Eindhoven, P. O. Box 513, 5600 MB Eindhoven, 
Nederland. 

1 B. Lorenz, R. L. Meng, and C. W. Chu, Phys. Rev. B 64, 
012507 (2001). 

2 V. G. Tissen, M. V. Nefedova, N. N. Kolesnikov, and 
M. P. Kulakov, Physica C 363, 194 (2001). 

3 K. Prassides, Y. Iwasa, T. Ito, D. H. Chi, K. Uehara, 
E. Nishibori, M. Takata, S. Sakata, Y. Ohishi, O. Shimo- 
mura, T. Muranaka, and J. Akimitsu, Phys. Rev. B 64, 
012509 (2001). 

4 E. Saito, T. Takenobo, T. Ito, Y. Iwasa, K. Prassides, and 
T. Arima, J. Phys.: Condens. Matter 13, L267 (2001). 

5 M. Monteverde, M. Nunez-Regueiro, N. Rogado, K. A. Re- 
gan, M. A. Hayward, T. He, S. L. Loureiro, and R. J. Cava, 
Science 292, 75 (2001). 

6 T. Tomita, J. J. Hamlin, J. S. Schilling, D. G. Hinks, and 
J. D. Jorgensen, Phys. Rev. B 64, 092505 (2001). 

7 T. Yildirim, O. Gulseren, J. W. Lynn, C. M. Brown, 
T. J. Udovic, H. Z. Qing, N. Rogado, K. A. Regan, 
M. A. Hayward, J. S. Slusky, T. He, M. K. Haas, P. Khal- 
ifah, K. Inumaru, and R. J. Cava, Phys. Rev. Lett. 87, 
037001 (2001). 

8 I. Hase and K. Yamaji, J. Phys. Soc. Jpn. 70, 2376 (2001). 

9 G. Gaspari and B. Gyorffy, Phys. Rev. Lett. 28, 801 (1972). 

10 K.-P. Bohnen, R. Heid, and B. Renker, Phys. Rev. Lett. 
86, 5771 (2001). 

11 D. Lampakis, A. Tatsi, E. Liarokapis, G. Varelogiannis, 
P. M. Oppeneer, M. Pissas, and T. Nishizaki, preprint 



cond-mat/0105447 



T. Muranaka, S. Margadonna, I. Maurin, K. Brigatti 
D. Colognesi, K. Prassides, Y. Iwasa, M. Arai, M. Takata 

and J. Akimitsu, J. Phys. Soc. Jpn., 70, 1480 (2001). 

13 J. D. Jorge nsen and D. G. Hinks, preprint cond- 



mat/0205486 . 
14 F. Marsiglio, Phys. Rev. Lett. 87, 247001 (2001). 
lo P. Ravindran, P. Vajeeston, R. Vidya, A. Kjekshus, and 

H. Fjelvag, Phys. Rev. B 64, 224509 (2001). 

16 P. Seneor, C.-T. Chen, N.-C. Yeh, R. P. Vasquez, L. D. Bell, 
C.U. Jung, Min-Seok Park, Heon-Jung Kim, W.N. Kang, 
and Sung-Ik Lee, Phys. Rev. B 65, 012505 (2002). 

17 S. Tsuda, T. Yokoya, T. Kiss, Y. Takano, K. Togano, 
H. Kito, H. Ihara, and S. Shin, Phys. Rev. Lett. 87, 177006 
(2001). 

18 P. Szabo, P. Samueli, J. Kacmarcik, Th. Klein, J. Mar- 
cus, D. Fruchart, S. Miraglia, C. Marcenat, and 
A. G. M. Jansen, Phys. Rev. Lett. 87, 137005 (2001). 

19 A. Junod, Y. Wang, F. Bouquet, and P. Toulemonde, to 
appear in: "Studies of High Temperature Superconductors", 
ed. A. Narlikar (N ova Science Publishers, NY.); preprint 



cond-mat/0106394 



On leave of absence from N. S. Kurnakov Institute of the 



A. Junod, private communication, October 2001. 

21 S. L. Bud'ko, G. Lapertot, C. Petrovic, C. E. Cunningam, 
N. Anderson, and P. C. Canfield, Phys. Rev. Lett. 86, 1877 
(2001). 

22 D. G. Hinks, H. Claus, and J. D. Jorgensen, Nature 411, 
457 (2001). 



5 



23 L. S. Schilling, J. D. Jorgensen, D. G. Hinks, S. Dee- 
myad, J. Hamlin, C. W. Looney, T. Tomita, to appear 
in: "Studies of High Temperature Superconductors ", ed 



A . Narlikar (N ova Science Publishers, N.Y.); preprint cond- 
m|at/0110268| 

24 A. Knigavko and F. Marsiglio, Phys. Rev. B 64, 172513-1 
(2001). 

25 J. J. Tu, G. L. Carr, V. Perebeinos, C. C. Homes, M. Stron- 
gin, P.B. Allen. W. N.Kang, E. M.Choi, H. J. Kim, and 
Sung-Ik Lee, Phys. Rev. Lett. 87, 277001 (2001). 

26 J. E. Hirsch, Phys. Lett. A 282, 392 (2001). 

27 J. E. Hirsch and F. Marsiglio, Phys. Rev. B 64, 144523 

(2001) . 

28 A. F. Goncharov, V. V. Struzhkin, E. Gregoryanz, J. Hu, 
R. J. Henley, Ho-kwang Mao, G. Lapertot, S. L. Bud'ko, 
and P. C. Canfield, Phys. Rev. B 64, R100509 (2001). 

29 H. Martinho, A. A. Martin, C. Rettori, O. F. de Lima, 
R. A. Ribeiro, M. A. Avi la, G. Pagliuso, N O. Moreno, 
and J. L. Sarrao, preprint cond-mat/0105204. 

30 P. Bordet, M. Mezouar, M. Nunez-Regueiro, M. Mon- 
teverde, M. D. Nunez-Regueiro, N. Rogado, K. A. Regan, 
M. A. Hayward, T. He, S. M. Loureiro, and R. J. Cava, 
Phys. Rev. B 64, 172502 (2001). 

31 V. A. Ivanov, M. van den Broek, and F. M. Peeters, Solid 
State Commun. 120, 53 (2001). 

32 V. A. Ivanov, J.J. Betouras, and F. M. Peeters, preprint 
cond-mat/0204085| . 

33 B. Lorenz, R. L. Meng, Y. Y. Xue, and C. W. Chu, Phys. 
Rev. B 64, 052513 (2001). 

34 I. I. Mazin and J. Kortus, Phys. Rev. B 65, R180510 

(2002) . 

35 T. Vogt, G. Schneider, J. A. Hriliac, G. Yang, and 
J. C. Abel, Phys. Rev. B 63, R220505 (2001). 

36 J. D. Jorgensen, D. G. Hinks, and S. Short, Phys. Rev. B 
63, 224522 (2001). 

37 J. Sak, Phys. Rev. B 6, 3981 (1972). 

38 E. Evans and D. L. Mills, Solid State Comm. 11, 1093 
(1972); Phys. Rev. B 8, 4004 (1973). 

39 Y. Kong, O.V. Dolgov, O. Jepsen, and O.K. Andersen, 
Phys. Rev. B 64, R0200501 (2001). 

40 V. A. Ivanov, Europhys. Lett. 52, 351 (2000). 

41 R. O. Zaitsev, JETP Lett. 73, 331 (2001). 

42 E. S. Choi and W. Kang, preprint :ond-mat/0104454 . 

43 K. P. Meletov, M. P. Kulakov, N. N. Kolesnikov, J. Arvan- 
itidis, and G. A. Kourouklis, JETP Lett. 75 406 (2002). 

44 K. Kobayashi and K. Yamamoto, J. Phys. Soc. Jpn. 70, 
1861 (2001). 

45 Feng-Yiang Li, Ru-Ju Wang, Shao-Chun Li, Liang- 
Chen Chen, Jia-Lin Zhu, Zhen-Xing Liu, Ri-Cheng Yu, and 
Chang-Qing Jin, Phys. Rev. B 65, 132517 (2002). 

46 V. P. Antropov, K. D. Belaschenko, M. van Schilfgaarde, 
and S. N. Rashkeev, to appear in: "Studies of High Tem- 
perature Superconductors" , ed. A. Narlikar (N ova Science 
Publishers, NY.); preprint |sond-mat/0107123 . 

47 M. Hanfland, H. Beister, and K. Syassen, Phys. Rev. B 39, 
12598 (1989). 



G 



